Calculations of spin-disorder resistivity from first principles 
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Spin-disorder resistivity of Fe and Ni is studied using the noncollinear density functional theory. 
The Landauer conductance is averaged over random disorder configurations and fitted to Ohm's 
law. The distribution function is approximated by the mean-field theory. The dependence of spin- 
disorder resistivity on magnetization in Fe is found to be in excellent agreement with the results 
for the isotropic s-d model. In the fully disordered state, spin-disorder resistivity for Fe is close 
to experiment, while for fee Ni it exceeds the experimental value by a factor of 2.3. This result 
indicates strong magnetic short-range order in Ni at the Curie temperature. 



The effects of thermal spin disorder on the transport 
properties of ferromagnets are of interest for both fun- 
damental and practical reasons. First, studies of spin 
disorder resistivity may provide important information 
on magnetic short-range order in ferromagnets, which 
is an intrinsic thermodynamic feature of spin disorder. 
Second, spin disorder intermixes the spin channels and 
thereby introduces finite spin diffusion length in ferro- 
magnets. Together with scattering on interface spin dis- 
order and spin-orbit scattering, this effect degrades the 
magnetoresistance of magnetic tunnel junctions and spin 
valves, which are used in magnetic field sensors and mag- 
netic random-access memories. 

A large amount of experimental data is available on the 
electric resistivity of ferromagnetic metals, most of which 
was collected many years agoi In these metals the resis- 
tivity has an "anomalous" contribution which contains 
signatures of magnetic phase transitions. This magnetic 
contribution /9 mag may be extracted from the measured 
resistivity assuming that the Matthiessen's rule is valid, 
and that the anomalous contribution is temperature- 
independent well above the Curie temperature^ 

The anomalous resistivity p ma g may be attributed to 
spin-disorder scattering, which was studied using the s- 
d model Hamiltonian by many authors This ap- 
proach assumes that 3d electrons in transition metals are 
localized at atomic sites, while the conductivity is due to 
mobile 4s electrons forming an itinerant band and cou- 
pled to the d electrons by exchange interaction. At finite 
temperatures the directions of the <i-electron spins fluc- 
tuate, and the conducting s electrons scatter from the 
inhomogeneous exchange potential. 

As the temperature is increased toward the Curie tem- 
perature T c , the spins become more disordered, and p m ag 
quickly increases, sometimes surpassing the phonon con- 
tribution. At elevated temperatures the s-d model cal- 
culations of Pmag were done using the mean-field approx- 
imation for magnetic thermodynamics. In this approx- 
imation, the spins are completely disordered above T c , 
and hence p m ag is constant. In the Born approximation 



below T c it declines as 

Pmag(T) = Pmag(T c ) [l - M 2 (T) /A/ 2 (0)] , (1) 

where M(T) is the magnetization at temperature T££& 
In general, p m ag is sensitive to magnetic short-range, 
rather than long-range order j! Quantitative comparison 
of the calculations of /O m ag(T) with experiment can thus 
provide valuable information on magnetic short-range or- 
der, which may be very strong in itinerant ferromagnets 
like Ni£ 

Unfortunately, since the s-d model assumes an 
isotropic conduction band, it can not predict the mag- 
nitude of Pmag(Tc) of specific materials. Moreover, treat- 
ing the 3d electrons as localized appears to be an inade- 
quate approximation for computing the spin-disorder re- 
sistivity, because the dominant contribution comes from 
interband (s-d), rather than intraband (s-s) scattering 
transitions. This conclusion was made by Goodings^ who 
solved the Boltzmann equation for the isotropic two-band 
model. Although this model is not applicable to transi- 
tion metals with complex Fermi surfaces, it clearly high- 
lights the sensitivity of the spin-disorder resistivity to the 
band structure. A parameter-free approach based on the 
density functional theory is clearly needed for quantita- 
tive comparison with experiment. 

In this paper we calculate p m ag for iron and nickel from 
first principles using the Kubo-Landauer formalism im- 
plemented in the tight-binding linear muffin-tin orbital 
(TB-LMTO) method within the atomic sphere approx- 
imation, which was extended to non-collinear systems. 
Our treatment of spin disorder is based on the adiabatic 
approximation, which assumes that the low-lying exci- 
tations of the magnet are well represented by the states 
where the magnetic moments at different atomic sites 
are rigidly rotated off of the magnetization axis»^ The 
state of the system is assumed to be the ground state 
constrained by these rotations. For a given noncollinear 
spin configuration, one can thus construct a TB-LMTO 
Hamiltonian or the Green's function matrix (which are 
no more diagonal in spin space) ^ We used the local 
spin density approximation and frozen atomic potentials 
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FIG. 1: The schematic picture of the system used in the cal- 
culations. Vertical lines indicate the embedding planes. 
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taken from the ferromagnetic ground state. The basis set 
included the states up to l max — 2 for both Fe and Ni (s, 
p, and d orbitals). 

The conductance was calculated using the principal- 
layer Green's function technique and the Kubo-Landauer 
formalism . 12 i 13 We consider a ferromagnetic layer FM(D) 
with spin disorder placed between two magnetically or- 
dered leads FM(O) made of the same material, as seen 
in Fig. [TJ As usual in the principal-layer technique, the 
"active region" consists of the FM(D) layer plus a few 
monolayers of FM(O) on each side. 

Because we use frozen atomic potentials to make cal- 
culations practical, care needs to be taken to ensure lo- 
cal charge neutrality. Indeed, one can view FM(D) and 
FM(O) as different materials which have different Fermi 
levels that must normally be matched by the contact 
potential. Frozen atomic potentials do not allow this 
to happen, and therefore we introduced a compensating 
constant potential shift in the disordered region which 
mimics the contact potential. This shift was found so 
that the charge in the central part of FM(D) averaged 
over disorder realizations is zero. Note that no mat- 
ter how the FM(0)/FM(D) interfaces are treated (self- 
consistcntly or not), they add contact resistances to the 
circuit. However, since the resistivity of the FM(D) ma- 
terial is extracted from the dependence of the resistance 
on the FM(D) thickness in the Ohmic limit, it does not 
depend on the interfaces, and therefore our crude treat- 
ment of the interfaces does not introduce any error. 

First we considered the paramagnetic state with no 
magnetic short-range order (completely uncorrclatcd 
spins with uniform angular distribution function). For 
the given thickness of the FM(D) layer, the resistance of 
the system was averaged over several disorder configu- 
rations (typically 12). As seen in Fig. [21 almost perfect 
Ohmic behavior was found for Fe down to the smallest 
thicknesses for the 4x4 lateral cell (16 atoms per layer). 
The convergence with respect to the lateral cell size was 
also checked. We found that for paramagnetic Fe /9 ma g 
is quite insensitive to this size, as seen in Table HI For 
Ni the thickness dependence of /0 mag is also fairly close 
to linear for the 2x2 lateral cell (8 atoms per layer). 

The calculated p mag for Fe (see Table HJ is about 20% 
greater than the experimental estimate.^ This result in- 
dicates that magnetic short-range order in Fe at T c is 
weak. However, for the similar model of fee Ni /9 mag 
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FIG. 2: The area-resistance product AR of the 
FM(0)/FM(D)/FM(0) systems as a function of the 
thickness of the FM(D) layer for bec Fe and fee Ni in the 
paramagnetic state. The lateral cell sizes were 4x4 for Fe 
and 2x2 for Ni, with edges along the [100] directions. The 
errorbars denote standard deviations of AR due to the finite 
number of spin disorder configurations. 



TABLE I: Spin-disorder resistivity p ma9 in /iQ-cm for para- 
magnetic bec Fe and fee Ni. The calculated values are given 
for different lateral cell sizes with edges along the [100] di- 
rections. The experimental values are the high-temperature 
asymptotes taken from Ref. 0- 
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is nearly 2.3 times greater than the experimental esti- 
mate. This large difference indicates either very strong 
magnetic short-range order, or a strongly reduced local 
moment at T ~ T c . 

Further, we calculated p m ag for ferromagnetic Fe in a 
range of temperatures. The spin configurations in the 
disordered Fe layer were generated using the mean-field 
theory, i.e. with an uncorrelated distribution function 
given by 



p(0) oc e - H ««-/V T t H cS (T) = 



3M(T)T C 
fj.M(0) 



(2) 



where 6 is the angle between the spin and the magnetiza- 
tion axis, y, is the magnetic moment at each site, M (T) 
is the magnetization at temperature T in the mean-field 
approximation, and H e g is the effective Weiss field. 

We used 4x4 lateral cells for all temperatures and 
thicknesses. The Ohmic limit (linear dependence of the 
conductance on length) was achieved for all temperatures 
down to about T c /3. This behavior agrees with what one 
would expect from simple mean-free path considerations. 
If we estimate the mean-free path as I = |Ai?bai/p, where 



3 



100 



80 



E 60 
o 

a 

<a. ? 40 



Calculations 
Experiment 



0.4 0.6 0.8 

A/(T)/A/(0) 



FIG. 3: Spin-disorder resistivity of Fe plotted against M 2 (T). 
The experimental data for p m a B (r) were taken from Ref. |2J, 
and those for M(T) from Ref. M- 



Ai?bai is the ballistic area-resistance product (this for- 
mula is exact for the free-electron model), we find that I 
does not exceed the lateral cell size in this temperature 
range. Therefore, quantum fluctuations of the conduc- 
tance are small, and the transport is diffusive. Another 
indication of the Ohmic behavior comes from the distri- 
bution of the current over the spin channels. The con- 
ductance is a sum of four partial conductances, Gr-ff , 
G11, G'it- Spin-conserving and spin-flip scattering have 
similar rates in our spin-disorder problem, and therefore 
the electrons lose memory of their spin over their mean- 
free path. Therefore, in the Ohmic limit the partial con- 
ductances must be proportional to the number of chan- 
nels in the left and right leads for the corresponding spin 
channels: G aa i oc M^M^,. This behavior was indeed 
observed down to T ~ T c /3. 

The results of p ma g(T) calculations for Fe are plot- 
ted in Fig. [3] versus M 2 (T) along with the experimental 
data taken from Ref. 0. Surprisingly, the calculated re- 
sults, which are based on first-principles electronic struc- 



ture and do not involve perturbation theory, show that 
the temperature dependence agrees with Eq. ([TJ) with a 
rather high accuracy. This striking agreement suggests 
that the temperature dependence of /3 mag is quite insen- 
sitive to the shape of the Fermi surface. 

In the low-temperature region (large M) the magnetic 
excitations are in reality dominated by spin waves, and 
the mean-field theory is inapplicable. Therefore, the fact 
that the theoretical and experimental curves in Fig. [3] 
have the same slope in this region is probably acciden- 
tal (the extraction of p ma g using Matthiessen's rule is 
also unreliable here). In the high-temperature region the 
comparison is more meaningful, and here we see that 
the calculated results are quite close to experiment at 
M 2 (T)/M 2 (0) ~ 0.5 and somewhat diverge as the tem- 
perature approaches T c . There may be several reasons 
for this behavior, including the dependence of the local 
moment on temperature which we neglected, and also 
the presence of magnetic short-range order. The resis- 
tivity is sensitive to short-range, rather than long-range 
order—, and the short-range order depends less strongly 
on the temperature. Therefore, the deviation of the the- 
oretical curve from experiment close to T c is expected. 
Quantitative analysis of these effects deserves further in- 
vestigation. 

In conclusion, we calculated the spin-disorder contribu- 
tion to resistivity for iron and nickel from first-principles. 
Comparison with experimental data suggests that mag- 
netic short-range order is weak in Fe and strong in Ni. 
The analysis of spin-disorder resistivity can thus pro- 
vide important information on magnetic short-range or- 
der in magnetic metals. We also studied the tempera- 
ture dependence of spin-disorder resistivity in Fe assum- 
ing mean-field thermodynamics. The resistivity depends 
linearly on M 2 (T)/M 2 (Q) in agreement with the predic- 
tions of the isotropic s-d model. 
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